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Abstract 

Bremsstrahlung of photons from highly relativistic electrons is inves- 
tigated. The cross section of the processes, which is suppressed due to a 
multiple scattering of an emitting electron in dense media (LPM effect) 
and due to photon interaction with electrons of a medium, is calculated 
with an accuracy up to " next to leading logarithm" and with the Coulomb 
corrections taken into account. Making allowances for a multiple scatter- 
ing and a polarization of a medium an analysis of radiation on a target 
boundary is carried out. The method of consideration of radiation in a 
thin target under influence of the LPM effect is developed. Interrelation 
with the recent experiment is discussed. 
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1 Introduction 



When a high-energy electron emits a soft photon via bremsstrahlung, the process 
occurs over a rather long distance, known as the formation length. If anything 
happens to an electron or a photon while traveling this distance, the emission 
can be disrupted. Landau and Pomeranchuk were the first who showed that 
if the formation length of bremsstrahlung becomes comparable to the distance 
over which a multiple scattering becomes important, the bremsstrahlung will 
be suppressed Migdal 0, developed a quantitative theory of this phe- 
nomenon. Side by side with the multiple scattering of emitting electron one has 
to take into account also an influence of a medium on radiated electromagnetic 
field. Since long distances are essential in the problem under consideration this 
can be done by introducing dielectric constant 6{uj). This effect leads also to 
suppression of the soft photon emission (Ter-Mikaelian effect, see in |^). A clear 
qualitative analysis of different mechanisms of suppression is presented in 
More simple derivation of the Migdal's results is given in |^. 

The next step in a quantitative theory of LPM effect was made in |^ . This 
theory is based on the quasiclassical operator method in QED developed by au- 
thors 0- One of the basic equations (obtained with use of kinetic equations 
describing a motion of electron in a medium in the presence of external field) is 
the Schrodinger equation in external field with imaginary potential (Eq.(3.3),[^). 
The same equation (without external field) was rederived recently in [|l^. The 
last derivation is based on the approach results of which is coincide basically with 
the operator quasiclassical method. In fill] a new calculation approach is devel- 
oped where multiple scattering is described with the path integral treatment. 

New activity with the theory of LPM effect is connected with a very successful 
series of experiments [|l2l - performed at SLAC during last years (see in this 



connection |T^). In these experiments the cross section of bremsstrahlung of soft 
photons with energy from 200 KeV to 500 MeV from electrons with energy 8 GeV 
and 25 GeV is measured with an accuracy of the order of a few percent. Both 
LPM and dielectric suppression is observed and investigated. These experiments 
are the challenge for the theory since in all the mentioned papers calculations 
are performed to logarithmic accuracy which is not enough for description of 
the new experiment. The contribution of the Coulomb corrections (at least for 
heavy elements) is larger then experimental errors and these corrections should 
be taken into account. 

In the present paper we calculated the cross section of bremsstrahlung process 
with term oc 1/L , where L is characteristic logarithm of the problem, and with 
the Coulomb corrections taken into account (Section 2 and Appendix A). This 
cross section is valid for very high energies when the LPM effect manifest itself 
for a photon energy of the order of an energy of the initial electron. In the 
photon energy region, where the LPM effect is "turned off', our cross section 
gives the exact Bethe-Heitler cross section (within power accuracy) with Coulomb 
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corrections. This important feature was absent in the previous calculations. The 
polarization of a medium is incorporated into this approach (Section 3). The 
considerable contribution into the soft part of the investigated in the experiment 
spectrum of radiation gives a photon emission on the boundaries of a target. 
We calculated this contribution taking into account the multiple scattering and 
polarization of a medium for the case when a target is much thicker than the 
formation length of the radiation (Section 4). In Section 5 we considered a case 
when a target is much thinner than the formation length. In this case the cross 
section has multiplicative form (probability of radiation times cross section of 
scattering for the given impact parameter) . In Section 6 a case of an intermediate 
thickness of a target (between cases of a thick and a thin target) is analyzed, 
polarization of a medium is not included. In Section 7 a qualitative picture 
of a spectral curve (an effective thickness of a target, position of a minimum) 
is discussed. In Section 8 we compare the theoretical curve for the intensity 
spectrum with the data. Although agreement between experiment and theory is 
rather satisfactory, an additional analysis should be done to obtain information 
about an accuracy of agreement between experimental data and theory. 

2 The LPM effect in an infinitely thick target 



As well known (see, e.g. |16|, Sec. 93) the formation length of radiation is (in this 



paper the system ^ = c = 1 is used) 

Iff' f 
lc=^, C = l + y^', e' = e-u;, 7 = -, (2.1) 

where e is the energy of the initial electron, uj is the energy of radiated photon, 
d is the angle between momenta of the photon and the initial electron. We 
consider first the case when the formation length is much shorter than thickness 
of a target l{lc ^ /). In this case the spectral distribution of the probability of 
radiation per unit time is given by expression (2.18), [0 (see also p|. Section 7.4) 



dW 

dui 



/CLT LkP' ( 

dTexp{-i—) ^-^¥?o(0,r) -i ( 1 + — 1 Vc^(0, 



(2.2) 



where a = = — — , functions ^^{^0, <p) satisfy an equation 

-Lo I 

^ - f A^.(x, r) = n(S(x) - S(0))^,(x, r) (2.3) 
with the initial conditions 

yPo(x, 0) = 5(x), (^(x, 0) = -«V5(x). (2.4) 
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Here n is the number density of atoms in the medium, x is the coordinate in 
two-dimensional space conjugated to the space (two-dimensional) of radiation 
angle t?, S(x) is the Fourier transform of the scattering cross section: 

S(x) =yc/¥exp(ixt?)a(^9), « = ^ = ^7- (2-5) 

For a screened Coulomb potential we have 

where "i^i = — , is the screening radius (a^ = O-SIobZ"^/^, - the Bohr ra- 

dius), K\ is the modified Bessel function. As we will show below, the main 
contribution to the probability is given by 

— ~'i9> — ^-(91 = — ^, where Ac = — = ( — ) is the electron Compton wave- 
X 7 m \mc J 

length. Expanding Ki{xdi) as a power series in xdi and introducing new vari- 
ables 

a [a 1 

i=o^' ^=VI^=~^' (2.7) 
z V 7 

we obtain for the spectral distribution of the probability of radiation 
dW 2a, 



du 7 



4 Re j dte-'' [i?i¥^o(0, t) + R2PV>{0, t)] , (2.8) 
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where Ri = — , R2 = — I , and the functions 09„ now satisfy an equation 

ee' e' e 

+ \n^ + 2C -l), Q = , C = 0.577... 

4 ^ m^uj 

with the initial conditions ipo{Q,0) = S{q), ifi{Q,0) = pS{q), the functions (po 
and ifi in ( |2.8| ) are rescaled according with the initial conditions (factors 1/7^ 
and 1/7'^, correspondingly). Note, that it is implied that in formulae ( p.2|) ,(|2.8| ) 
subtraction at = is made. 

The potential V{q) ( |2.9| ) corresponds to consideration of scattering in the 
Born approximation. The difference of exact as a function of Za potential V{q) 
and taken in the Born approximation is computed in Appendix A. The potential 
V{q) with the Coulomb corrections taken into account is 



V{q) = -QQ^[ln^^'dl + ln^ + 2C -l + 2fiZa 



(2.10) 
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where ■i?2 = exp(/ — 1/2), the function / = f{Za) see in ( [A. 101 ). 

In above formulae Q is space of the impact parameters measured in the Comp- 
ton wavelengths Ac, which is conjugate to space of the transverse momentum 
transfers measured in the electron mass m. An operator form of a solution of 
Eq. dU) is 

MQ^t) = exp{-iHt)^o{Q,0) =< Q\exp{-iHt)\0 >, H = -iV{Q), 
(p{g,t) = exp{—iHt)p(po{g,0) =< ^| exp(— ii/t)p|0 >, 

(2.11) 

where we introduce the Dirac state vectors: | ^ > is the state vector of coordinate 
Q, < q\0 >= S{q). Substituting ( p.ll|) into (p.8| ) and taking integral over t we 
obtain for the spectral distribution of the probability of radiation 

^ = ^Im < 0|i?i {G-' - Go + i?2P - Go p|0 >, (2.12) 

where 

G = p2 + l-iV, Go = p^ + l. (2.13) 

Here and below we consider an expression < 0|...|0 > as a limit: lim x ^ 0, 
lim x' — > of < x|...|x' >. 

Now we estimate effective impact parameters Qc which give the main contri- 
bution into radiation probability. Since characteristic values of Qc will be found 
straightforwardly at calculation of ( p.l2|) , we estimate characteristic angles -^c 
connected with Qc by an equahty = 1/ ili^c) ■ The mean square scattering 
angle of a particle on the formation length of a photon Ic ( p.l|) has the form 



,5.1^^„U„ « ^i|.„ « . (2.14) 

When -t^l <^ 1/7^ ^he contribution in the probability of radiation gives a region 
where ( ~ l{{}c = 1/7)5 ^^is case Qc = I. When {}s ^ I/7 the characteristic 
angle of radiation is determined by self-consistency arguments: 

- € - % = 4^1n4^' 7?1^4l2=l' 4Qg^ln ^ = 1. (2.15) 

7 Cc7 Y^i Cc 7^1 Y^iQi 

It should be noted that when characteristic impact parameter becomes smaller 
than a radius of nucleus the potential V^(^) acquires an oscillator form (see 
Appendix B, Eq.O) 

Vie) =Qg'(^ln^ - 0.208^ (2.16) 

Allowing for estimates (|2.15|) we present the potential V{e) (|2.9|) in the fol- 
lowing form 

Vie) = Vcie) + v{e), K(e) = q = QL, l = in \ ^ , 

2 / .2 X ^ "^'^^ (2.17) 



^(e) = -^(2G + ln.^ „ 
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The inclusion of the Coulomb corrections f{Za) and -1 into Iwdl diminishes 
effectively the correction v{q) to the potential Vc{q). In accordance with such 



division of the potential we present propagators in expression ( p.l2| ) as 

- = G-^ - G-^ + - G^^ (2.18) 

where 

G^ = p'^ + l-tV„ G = p2 + 1 - iV^ - IV 

This representation of the propagator G~^ permits one to expand it over "per- 
turbation" V. Indeed, with an increase of q the relative value of the perturbation 

V 1 

is diminished (— ~ — ) since effective impact parameters diminish and, corre- 



spondingly, the value of logarithm L in ( p.l7| ) increases. The maximal value of 
L is determined by a size of a nucleus Rn 



^~21n^ 

m A2 



Lma. = In ^ ~ 2 In ^ = 2Li, (2.19) 



where 0^2 = asexp(— / + 1/2). So, one can to redefine the parameters and {}i 
to include the Coulomb corrections. 

The matrix elements of the operator G^^ could be calculated explicitly. The 
exponential parametrization of the propagator is 

00 

G-^ = ij dte-'^exp{-iH^t), H^ = p^- iqg^ (2.20) 


Below we will use matrix elements of the operator exp{—iHct) 

< Q^\exp{-iH^t)\Q^ >= K^{Q-^,Q^,t). (2.21) 



The function Kc{qi, Q2,t) satisfies the Schrodinger equation ( p.9| ) over each of 
two (symmetrical) variables Qi and Q2 with V = qg^ and the initial condition 

i^c(ei, ^2,0) =5(^2-^1)- (2.22) 
We will seek a solution in the form (see also |^) 

Kc{Qi,Q2,t) =exp[a{t){Ql + qI) + 2p{t)Q^Q2 + j{t)_ . 
Substituting this expression into ( pl9|) we find a set of equations for a,/?, 7 

a = Aia^ - q, /? = 4ia(3, 7 = 4ia. (2.23) 



The initial conditions for this set follows from definition ( 2.21 
lim < Qi\exp{—iHct)\g2 >^< Qi\exp{—iHQt)\Q2 >= 



t-^ 



^ /rfVexp(z(^2-^i)p-^P^t)=^exp[^i^i^^) (2.24) 



(2vr)2 

= Ko{Q2, Qi,t), 
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where Hq = p^. From ( |2.24| ) one has the initial conditions at t — > 

-ln(47rit). 

The solution of the set (|2.23|) satisfying these initial conditions is 



^^^^ 



= — cothi/t, (3{t) 



7(t) = -ln(sinh.t)+ln£., 



(2.25) 



(2.26) 



where v = 2\Jiq. As a result, we obtain the following expression for the sought 
function 



V 



Alii sinh ut 



[q\ + qI) cothz/t 



sinh vt 



Q1Q2 



}. (2.27) 



Substituting formulae ( |2.20| ) and ( p. 271 ) in the expression for the spectral distri- 
bution of the probability of radiation ( p.l2|) we have 



a 



27r7" 



rim $(//) 



$(i/) = u dte 



Ri 



f 1 



Vsinh 



z z 



iuRo 



sinh z 



(2.28) 



where z = ut. This formula gives the spectral distribution of the probability of 
radiation derived by Migdal 0. However, here Coulomb corrections are included 
into parameter v in contrast to 0. 

We now expand the expression G^^ — G^^ over powers of v 



G 



-1 



G 



-1 



G-\-w)G-' + G-\~tv)G-\-tv)G-' + ... 



(2.29) 



Substituting this expansion in ( |2.18| ) and then in ( |2.12| ) we obtain decomposi- 
tion of the probability of radiation. Let us note that for Q 1 the sum of the 

probability of radiation ( 2.28|) and the first term of the expansion ( |2.29|) 
gives the Bethe-Heitler spectrum of radiation, see below (|2.39|) . At Q > 1 the 

expansion ( p.29|) is a series over powers of — . It is important that variation of 

Ij 

the parameter by a factor order of 1 has an influence on the dropped terms in 
only. 

In accordance with (|2.18 ) and ( 2.29 ) we present the probability of radiation 
in the form 

dW dWc dWi dW2 

+ ^^ + ^^ + ... (2.30) 



dbJ 



duj 



duj 



du 
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The probability of radiation — — ^ is defined by Eq.( |2.28D . In formula (|2.12| ) with 

dui 

allowance for (|2.18| ) there is expression 

oo oo 

-i < 0|G"^ -G;^|0 >= J dtij dt2e-'^'''+'^^ J d^gK^{0, Q,ti)v{Q)K,{Q,0,t2) 



oo oo oo 

+ JdhJ dt2 J dt3e-'^''+'''+''^ J d^Qi J d^02Kc{O, Q„h)v{Q,)K,{Q„ Q^,t2) 





Xv{Q2)K^{Q2,0,t3) + ... 



dWi . 



(2.31) 



where the matrix element Kc is defined by ( p.27| ) . The term in (|2.30|) 

corresponds to the first term (linear in v) in (p.31|) . Substituting ( 2.27 ) we have 



dWi _ 2a 
duj 7 



oo oo 



^Re J dhj dt2e~"^''+'''> J d^gvio) 



X exp 







qg , _ , . 
(coth uti + coth 1/12) 



Ri + 



vr^z/^ sinh uti sinh z/t2 
Aq^g^ 



sinh z/ti sinh z/t2 



Ri 



(2.32) 



where v = 2^iq. Substituting in ( p.32|) the exphcit expression for vi^Q) and 
integrating over di^g and d{t\ — ^2) we obtain the following formula for the first 
correction to the probability of radiation 



dWx 
duj 



a 



Im F{v)- F{v) = [ [RMz) - 2^i?2/2(^)], 

J smh z 



fi (z) = ( In + In — — In sinh z — C ) g{z) — 2 cosh zG{z) 



sinh z 



, (7(2;) = z cosh z — sinh 2;, 



(2.33) 



G{z) = y (1 — ycothy) dy, t = ti + t2, z = ut 



As it was said above (see ( |2.15| ), ( |2.19| )), gc = 1 a.t 



l{q = QLi). If the parameter z/i > 1, the value of gc is defined from the equation 
( |2.15|) , where "i^i ^^2, up to gc = Rn/Xc- Then one has 

\ngl + \n'^ = hn{gt4QL)-z^ = -z^, g^AQL = 1. (2.34) 

It follows from ( p.34| ) that expression (|2.15|) for which we chose a priori, 
corresponds to the mean value of g'^. From the above analysis we have that the 
factor at g{z) in expression for fi{z) in ( p.33| ) can be written in the form 

{Ingl + In — — Insinhz — C) — * (lni^o'*^(l — 1^0) ~ ~ Insinh^; — C), (2.35) 



where 



Z/q = \v\ 



Aq = AQL{q^ 



(2.36) 



is the Heaviside step function. 
When a scattering is weak (z/i ^ 1), the main contribution in (|2.33| ) gives a 
region where 2; -C 1. Then 



A(;.)^-(C + ln(^t))- + -z 



3 ^ ,2 



3 '3 



C-ln(2t)), 



-Im F{u) = -Im {R2 - Ri) , L Li. 
The corresponding asymptotes of the function ^{v) (|2.28|) is 



(2.37) 



$(z/) 



"6 



(i?i + 2i?2), (kl<l) 



(2.38) 



Combining the results obtained ( |2.37|) and ( p.38|) we obtain the spectral distribu- 
tion of the probability of radiation in the case when scattering is weak (|z/| ^ 1) 



dW dW, 



duj duj 
_ a 2Q 

~ 27172 T" 



+ 



a 



duj 
Ri (Li 



27172 
1 



Im 



2R2 Li 



3m2 



UJ 



ln(l83Z-i/=^) -^-/(Za 



(2.39) 



+2 (^1 + '-Pj (in (183Z-V3) + 1 _ /(z«)) , 

where Li is defined in ( p.l9|) . This expression coincide with the known Bethe- 
Heitler formula for probability of bremsstrahlung from high-energy electrons in 
the case of complete screening (if one neglects the contribution of atomic elec- 
trons) written down within power accuracy (omitted terms are of the order of 

powers of — ) with the Coulomb corrections, see e.g. Eq. (18.30) in 0, or Eq.(3.83) 



m 
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The integral in the function Im F(z^) ( p.33| ) which defines the first correction 
to the probability of radiation (|2.33|) can be transformed into the another form 
containing the real functions only 



-lmF{u) = Di{uo)Ri + -D2iuo)R2; s = -=-, 

s \J2vq 

00 



^i(^^o) 



dze 



sinh z 



d{z) sin sz + -^9{z) cos sz 



X 



1 



diz) - -g{z) 



(sin sz + cos sz) + -^g{z) (cos sz — sin sz) > 



dze 



sinh z 



(2.40) 



d{z) = (In I'o'&i^ — ^0) — In sinh z — C)g{z) — 2 cosh zG{z) 
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where the functions g{z) and G{z) are defined in ( |2.33| ). The form ( |2.40D is 
convenient for numerical calculations. Note, that parameter s in ( p.40|) is two 
times larger than used by Migdal [0. 

At z/q ^ 1 the function F^u) (see ( p.33| ) and ( |2.35| )) can be written in the 
form 



Fiu)- 

Integrating over z we obtain 



-^[RMz)-2tRMz)]. 
smh z 



-Im F{u) = lR, + ^i^\n2-C + ^^ R,. 

Under the same conditions {uq ^ 1) the function Im $(i^) ( ^.281 ) is 

Im 



-R^ + -j^R,. 



Thus, at z/q ^ 1 the relative contribution of the first correction 
by 



dWi 



dWi _ 1 
dW, ~ 2L(^ 



TT 



ln2 - C + - ) ~ 



0.451 



(2.41) 



(2.42) 



(2.43) 



is defined 



(2.44) 



where L{gc) = In 

In the above analysis we did not consider an inelastic scattering of a projectile 
on atomic electrons. The potential Ve{Q) connected with this process can be 
found from formula ( |2.10| ) by substitution Z"^ Z.'di ^ 'de = 0.153^91 (an 
analysis of an inelastic scattering on atomic electrons as well as the parameter 
•de can be found in JT^)- The summary potential including both an elastic and 
an inelastic scattering is 



1 



V{q) + V,{q) = -Q(l + -)Q'[\nY^', + In ^ + 2C + 

Q 



Z+l 



(2.45) 



where 



Qef = Q{1 + — ), ^ef = A exp 



iZfiaZ) - 1.88) - - 



3 An influence of the polarization of a medium 

When one considers bremsstrahlung of enough soft photons lo < ujol, one has to 
take into account the effect of a polarization of the medium. In a dense medium 
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the velocity of a photon propagation differs from the hght velocity in the vacuum 
since the index of refraction n{uj) 7^ 1 

Because of this the formation length diminishes as well as the probability of 
radiation (see Q], the qualitative discussion may be found in 0). For analysis 
we use the general expression for the probability of radiation, see Eq.(2.1), |Q. 
The factor in front of exponent in this expression (see Eq.(2.2), P]) contains two 
terms A and B, the term A is not changed and the term B contains combination 

V ^t? + n-^, Ko = , 3.2 

UO 27^ UJ 

and its dependence on (term of the order 1/7^) may be neglected also. With 
regard for the polarization of a medium the formation length ( p.l|) acquires a 
form 

2" 



UO 



27^ ' ^ 



1 + 7^^^+^^"^°^ 



. (3.3) 

\ UJ J 

So, the dependence on uq manifests itself in the exponent of Eq.(2.1), and 
respectively in the exponent of ( |2.2| ) only: 

a 2 (1 — —v] ~ OK = a, k = 1 + k^. (3.4) 

e — UJ \ UJ J 

Performing the substitution a — > a in formula (p.7|) we obtain for the potential 
Vc{,~q) = qg"^, q = QL{qc), Q = ^, L{Qc) = In 2^2-2 ^ 

(3.5) 

The substitution ( ^^) in the expression for the probability of radiation (|2.8|) 
gives 

R\ — > i?!, R2 — -^2^^ = R2- (3.6) 

The value of the parameter Qc in (|3.5|) is determined by equation (compare with 
Eq.O)) 

Aq,^QL{q,) = 1, for 4QL(1)>1. (3.7) 

In the opposite case Qc = 1 and this is possible in two intervals of the photon 
energy uj\ 
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1. for <^ 1 when the multiple scattering and effects of the polarization of 
a medium are weak; 

2. for kq ^ 1 when effects of the polarization of a medium become stronger 
then effects of the multiple scattering (z/q < n). 

In an intermediate region we substitute q^k in Eq.( p.7|) . After it we obtain 

the equation for which coincides with Eq. ( p.34| ), see also ( |2.36| ): 



ire 



(3.8) 



Thus, for <1 we have 



while for z/g < 1 we have 



Qc^ 



qIk 



K 



(3.9) 



^0 



4gL(i) 



K 



gin 



(3.10) 



The spectral distribution of the probability of radiation ( p.39| ) with allowance for 
polarization of a medium have the form 



dW 



a 



du! ^ny 



-Im 



2L(^c) 



(3.11) 



where 

<l(i?l, i?2) = $ (i?l, i?2), F(i?i, i?2) = F(i?l, i?2), 

We consider now the case when an influence the polarization of a medium man- 
ifests itself in the conditions of the strong LPM effect (t'o ^ !)• This influence 
becomes essential for low energy photons, when the mean square angle of the 
multiple scattering (|2.15|) on the formation length of a photon becomes smaller 

than ujI/oj"^ [vq = — < 1, Kq ^ 1 ). Indeed, in the case z>o ^ l(z^o ^ ^^o) 



can use asymptotes of functions $(z^) and F^p) at z/q ^ 1 (see ( ^.41| ), ( |2.43| )), 
we have 



dW a ^ + 

— — = -K2H 

duj 

r = ~ . . , L{g, 



V2 



= ln."- 



(3.12) 



\2^2 ' 



In the opposite case <^ Kq, the characteristic momentum transfer in the used 
units {(c) are defined by value K^g^ = 1), one can use asymptotic expansions 
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( |2.37D and ( p^.38D and we have for the spectral distribution of the probability of 
radiation 

^p + TT;- fyZa) = — Lp + — - f{Za^ 



dcu 3 m'^ounl V ^ 12 ' ^ V 3n m-f^ V ^ 12 

(3.13) 

where f{Za) is defined in ( |A.1(J| ), Lp = In ^183Z~^'''^ko^ . The results obtained 
agree with given in where calculations are fulfilled within logarithmic accu- 
racy and without Coulomb corrections. It is seen that a dependence of spectral 
distribution on photon energy [udu) differs essentially from the Bethe-Heitler 
one {duj/uj), the probability is independent on density n. 

The formula (|3.13| ) is applicable only up to value kq = Xc/ Rn or if > Ub, 
where 

u;, = ^ujol ^ aZ^/'luJo; - > aZ^"'-. (3.14) 
Ac em 

For example, for electrons with energy £ = 25 GeV and gold target (c^o = 80 eV^ 
one has lo^, ~ 125 KeV . For lo < lo^ one has take into account the form factor of 
a nucleus (see Appendix B). In this case the argument of the logarithm in (3.14) 



ceases its dependence on photon energy oj. In the limit uj <^ uji, the spectral 
distribution of the probability of radiation is 



dW 4 Z'^a^uj 



In -^-0.02 (3.15) 



dijj 37r 77172 V 

4 A target of a finite thickness 

In the case when a finiteness of a target is essential the probability of radiation 
is defined not only by the relative time r = ^2 ~ ^1 as in Section 2. The used 
radiation theory is formulated in terms of two times (see eqs.(2.1) - (2.3) of p). 
Proceeding from this formulation we can obtain more general expression which 
takes into account boundary effects. With allowance for polarization of a medium 
we have for the spectral distribution of the probability of radiation 

, , 00 t2 

dw 4q; t t 

— = — Re y dt2 J dtiexp{-ifi{t2)t2 + ifi{ti)ti) ^^^^ 

—00 —00 
X [riV9o(0, t2, h) - ir2Vcp{0, t2, ti)]. 



where 



In liiiTn'^ 

/i(t) = 7?(-t) + 79(r -t) + nm^iT -t), T=- = -— 
ri = —, r2 = 1 + k = 1 + Kq, 



(4.2) 



here / is the thickness of a target, kq is defined in (|3.2|). So, we split time interval 
(in the used units) into three parts: before target (t < 0), after target (t > T) 
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and inside target (0 < t < T). The functions ip^{Q,t2,ti), ip^j, = (p^{(po,(p) 
satisfy the equation ( |2.9| ), but now the potential V depends on time 



<^o(e, ^1, ^i) = S{q), (p{q, tl, tl) = pS{q). 

Using an operator form of a solution of Eq. ( |4.3| ) (compare with ( p.ll| )) we can 
present the probability ( [4.1|) in the form 

t2 



-Re / dt2 / dtiexp {—ifi{t2)t2 + ifi{ti)ti) 

-i J n{t)dt 



duo UJ 

-OO -OO ^ ^ (44) 

X (0|ri5(t2, tl) + r2p5(t2, ti)p|0) , 5(t2, tl) = T exp ' 



where the symbol T means the chronological product. Note, that in ( [4.1D and 
([4.4|) it is implied that subtraction is made at \^ = 0, /i(t) = 1 (k = 1). 
Integrals over time in (4^) we present as integrals over four domains: 

1. tl < 0, < t2 < T; 

2. < tl < T, < t2 < T; 

3. < tl < T, ta > T; 

4. tl < 0,t2 > T; 

in two more domains ti.2 < and ti.2 > T an electron is moving entirely free and 
there is no radiation. We consider in this Section the case, when the thickness 
of a target L is much larger than formation length ( p. 3D or (z/q + n)T S> 1. In 
this case domain 4) doesn't contribute. The contributions of other domains are 

OO ^ ^ 

h- dti / (it2exp (i(ti - Kt2)) exp {-iHt2) exp (iHoti) = -— — — 7^——, 

-00 

T t2 00 

/2 = y rft2 J (itiexp + K)(t2 - tl)) ~ T y rfr exp + K)r) 





T 1,11 



/ Tfirexp {—i{H + k)t) = ~i— h 

J H + K 



H + n {H + nf " Hq + IH + k 

(4.5) 

where Hq = p^. The term in I2: —iT/{H + k) describes the probability of 
radiation considered in previous Sections. All other terms define the probability 
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of radiation of boundary photons [|. So, making mentioned subtraction we have 
for the spectral distribution of the probabihty of radiation of boundary photons 

^ = —Re (0|riM + r2pMp|0) , M = M^/^ + M^^^ + Mq; 

^ 2 2 1 1 



Mn 



(1) 



(p2 + /t)(p2 + l) (h + k){p^ + 1)' ^ {H + K 
1 2 1 



[p- 



(p2 + ft:)(p2 + l) ^ (p2 + fi:)2^ 



(4.6) 

For a convenience here we made the subtraction in two stages: first 
(in My\ My^) we subtracted terms with V = and second (in Mq) we sub- 
tracted terms with both V = and k = 1. 

We consider important case when both the LPM effect and the polarization 
of a medium are essential. We will calculate the main term with V{q) = Vc{q), 
see ( p.l7|) . Needed combinations are 











oo oo 

- J dti J dt2exp{-i{ti + Kt2)) J (fgKc{0,Q, t2)Ko{Q,0,ti] 









1 



{H + ny 
(0|pMop|0) = 



tdtexp{-iKt)KciO,0,t), 





oo 



(2vr) 



1 + 



iMolO) = -, 



In K 



(4.7) 



where the functions Kq{q2, Qi,t) and Kc{Q2,Qi,'t) are defined in ( |2.24|) and 
(|2.27|) . Substituting into (|4.7| ) the explicit expressions for these functions, calcu- 
lating the vector derivatives as indicated in ( |4.6| ) we have for contribution of the 
first term in 



(1) 



duj 



2a 
■noj 



oo oo 

r2 Re iJ^ j dti j dt2 exp {—i{ti + ^^2)) 



1 



(z/ti + Z/t2 

2a 



—r2 Im u [ dti [ dt2 exp {-i{ti + 12)) 

TILU J J 



1 



(sinh z/t2 + ^^1 cosh z/t2 



1 



tanh z/t2 + i^ti vt2 + vt\ 



(4.8) 

where i) = v/k, the second term in the square brackets is the subtraction term 
in accordance with ([4.6|) (the term My"*). For practical use it is convenient to 

^Radiation of boundary photons in an inhomogeneous electromagnetic field was considered 

in 
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write the probability 
can be written as 



using real variables. After some transformations it 



dw. 



(1) 



where s 



_ 2a 
1 

t-y + y/K 
1 



r2 J dtexp{—t) (cost + sin t) J dy 



t — y + s tanh {y / ks) 



(4.9) 



V2uo 

tions with the second term in 
L6|)) we have 



parameter z/q is defined in ( p. 36 ). Repeating the same opera- 
(this is the contribution of the term My in 



dw, 



(2) 



du 
a 

7TUJ 

a 



= -^r2 Re / tdtexp {—iKt) 



oo 

r2 Re J dzexp (^—i—^ 



oo 

r2 J dzexp {—sz) cos sz 



sinh 



sinh 2; z 



(4.10) 



sinh 2 z 



where z = ut, s 



is calculated in ( [4.7[ ) 



i^o = '^o/'^- The contribution of the term Mq in (^.6|) 



(3) 



duj 



a 

TTUJ 



1 + 



lnK-2 



(4.11) 



The complete expression for the spectral distribution of the probability of ra- 
diation of boundary photons, in the case when both the LPM effect and the 
polarization of a medium are taken into account, is 



dwb 
doj 



3 rfu;f 



E 

k=l 



duj 



(4.12) 



We consider now the limiting case when LPM effect is very strong (z/q ^ 1). 
In this case we find for probabilities in formulae ( [4. 9]) and ( |4.1CI|) 



dw, 



(1) 



dw 



dw 



(2) 



duj 



2a 

71UJ 

a 

TIUJ 



1 ~ ^ TT 

In z^o - 7 + 



n - 1 8v^i>o y2i>o 
2 



TT 



In i>o + 1 - C + - 



1 - In 2z>o + C 



TT 



6y2i>o, 



(4.13) 
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Substituting asymptotes obtained and ( |4.11| ) into ( [4.12D we have 



dwb a \ 

-r- = — <ri+r2 

aUJ TTUJ 



Inz/o - 1 - C -ln2 



V2( 

K 

1^0 \ 24 



+ In I/O + 1 



(4.14) 

As one can expect, the probabihty of radiation at z/q ^ 1 + /^o depends on the 
polarization of a medium in the term oc 1 /z/q only. 



In the opposite case z>o ^ l(z/o ^ 1), the probabilities dw^^\ dw)^' oc Vq and 
probability of radiation of boundary photons is determined by the polarization 
of a medium. Just in this case radiation of boundary photons is known as the 
transition radiation: 



,(2) 



dwh dw 



(3) 



duj 



duo 



a 

TfUJ 



ri + r2 



In K - 2 



In the case of weak LPM effect z/q ^ 1 (cj ^ e) we have 



dwb 
duj 



a 



-r2 
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(4.15) 



(4.16) 



In this case what we calculated as the boundary photons contribution is actually 

correction (very small) to the probability I— — ( |2.39| ) which in this case has 

16 ^ 
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which follows from the decomposition 



additional (suppression) factor 1 

of the function Im $. 

The LPM effect for the case of structured targets (with many boundaries) 

The radiation of the boundary photons with regard 
2^] (for uj <ti e), the polarization 



was analyzed recently in 
for the multiple scattering was considered in 



of a medium was added in |2l|] and [^. Our results, which are consistent with 
obtained are presented in more convenient for application form and the 

Coulomb corrections are included. In these papers the probability of radiation of 
boundary photons (under condition of applicability of Eq. (^4.14 )) was analyzed 
also to within the logarithmic accuracy (see Eq.(20) in |^ and Eq.(15 in p2|)). 
This accuracy is insufficient for parameters connected with experiment |ll2|-|l4 . 
For example, for e = 25 GeV and heavy elements the value z/q equates k for 
z/q ~ 20. One can see from Eq. ([4.14|) that in this case In z/q is nearly completely 
compensated by constant terms. 



5 A thin target 

Finally we consider a situation when the formation length of radiation is much 
larger than the thickness / of a target (a thin target, I). In this case the ra- 
diated photon is propagating in the vacuum and one can neglect the polarization 
of a medium. 
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Operator S{ti,t2) ([4. 4)) we present in the form 



S{t2,ti) = Texp 



t2 

-i J n{t)dt 



exp {-iHot2) C{t2, h) exp {iHoti 



(5.1) 



C{t2,ti) = exp {iHot2) S{t2,ti)exp {-iH^ti). 
Differentiating the operator C{t2,ti) over the first of arguments we obtain 



dt 



exp (iHot) V{q, t)S{t, ti) exp {-iH^ti) = -V{q + 2pt, t)C{t, t 



1), 



(5.2) 

where V{Q,t) = V{Q)g{t) (see ( [^.9D , ( |4.3| )). The formal solution of this equation 
with the initial condition C{ti,ti) = 1 has the form 



^ih^ti) = Texp 



t2 

J dtV{Q + 2pt,t) 



(5.3) 



where T means the chronological product. This solution is exact. Now we take 
into account that we are considering a short characteristic time contributing into 
integral ( p.3|) , or more precisely 



^ , , 2^1 

t<T = -, l<^l,= — T<-, 

2 aC C 



(5.4) 



where Ic-, C ^ire defined in ( p.l|) . Since the main contribution give p ~ a/C, 

^ ~ pt l/\/C ~ Q , where p is characteristic mean value of operator |p|, 

one can neglect by the term 2pt in (pl3|), so that 



C{t2,ti) ^ exp 



i2 



dtV{Q, t) 



ti 



(5.5) 



In the probability of radiation enters the expression (cp (p.l2|) , ( [4.4|) ) 
(0 |exp i-iHot2) {C - 1) exp {iH^h)\Q) 

= J d^giC-l){0\exp{-tHot2)\Q){Q\exp{tHoh)\0). (^-6) 



Using an explicit form ( [2.24] ) of the matrix element (0 |exp (— i/fo^2)| q) and ne- 



glecting terms of the order ~ T{l/lc) one obtains starting from (f4.4|) for the 
spectral distribution of the probability of radiation 



dw 



th 



a 



dt 



du An'^LU 



X exp 

a 



(exp(-yT) - 1) 



d'g \riK',{g) + r2K',{g)] (1 - exp{-VT)), 



(5.7) 
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where pi (P2) is the operator p = — acting on the function of /ti {0^/12), 
Kn is the modified Bessel function. Here we took into account that in our case 
contribute domain \t2\ ^ T and ti < 0,^2 > since in domains ti^2 < and 
ti 2 > T an electron is moving entirely free and there is no radiation. In implicit 
form the factorization contained in (|5.71 ) is presented in [23|. If 
V{q = 1)T <C 1 one can expand the exponent (the contribution of the region 
^ 1 is exponentially damped because in this region Kq i^q) oc exp{—g)). In 
the first order over VT using the explicit expression for the potential ( p. 101 ) we 
have to calculate following integrals: 



(5.8) 



00 1 °° 1 1 

/ K^oi0)0'dg = -, J KI{q) In qq^cLq = -{\n2-C+- 



00 00 

j Kl{Q)Q^de = 3, / KI{q) In QQ^dQ = -{\n2-C-- 



Substituting these integrals one obtains in this case the Bethe-Heitler formula 
with the Coulomb corrections ( p.39| ). 

We analyze now the opposite case when the multiple scattering of a particle 
traversing a target is strong {y{Q = 1)T ^ 1, the mean square of multiple 
scattering angle > 1/7^)- We present the function V{g)T (see (^), (|2.10|) 
and ( CT )) as 



V{o)T = L^-2C]= A,^ In ^ = A,^ f In ^ - In ^ 




2 / ' 



(5.9) 

where gt is the lower boundary of values contributing into the integral over g. 
Substituting this expression into ( |^ ) we have the integral 



00 f 
2% j gdgKl{g) j 1 ~ ^xp 



-kg^ [ 1 - ^In ^ 
V Lt gf^ 



ttJ. 



(5.10) 



In this integral we expand the exponent in the integrand over 1/Lj keeping the 
first term of the expansion. We find 



J = Ji + J2, Ji = 2 J K'l{g) 1 — exp y—kg^j gdg 



00 

= 2k j dgg"^ [KQ{g)K2{g) - Kl{g)\ exp (-kg"), 


2k y 

J2 = —y- / K'l{g) exp {~kg^\ In —g^dg 
Lt J ^ ■' Pt 



(5.11) 
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In the integral Ji we performed an integration by parts. In the integrals in ( |5.11| ) 
it is convenient to substitute z = kg^ then 





1 

J2 = -—jKl[^^eM-z)\^zzdz. 



exp(— 

(5.12) 



Expanding the modified Bessel functions Kn{x) at x <^ 1 and taking the integrals 
in the last expression we have 



C_ 

t 



k = "^^nl {Lt + l- 2C). 

'2 



(5.13) 



In the term with Kq in (|5.7] ) the region ^ ~ 1 contributes. So we have 

oo oo 

Js = 2! Klio) (1 - exp(-\/T)) QdQ-2j K^{g)gdg = 1. (5.14) 



Substituting found J and J3 into (|5.7|) we obtain for the spectral distribution of 
the probability of radiation in a thin target at conditions of the strong multiple 
scattering 

— — = — (n+raJ). (5.15) 

The logarithmic term in this formula is well known in theory of the collinear 
photons radiation at scattering of a radiating particle on angle much larger than 
characteristic angles of radiation ~ I/7. It is described with logarithmic accuracy 
in a quasi-real electron approximation (see p4|, Appendix B2). 



The formula ( p.7| ) presents the probability of radiation in the case when the 
formation length Ic ^ /. It is known, see e.g. 0, that in this case a process of 
scattering of a particle is independent of a radiation process and a differential 
probability of radiation at scattering with the momentum transfer q can be 
presented in the form 

dW^ = dws{q)dwr{q,^)^ (5.16) 

where dws{ci) is the differential probability of scattering with the momentum 
transfer q which depends on properties of a target. The function dwr{ci,k.) is 
the probability of radiation of a photon with a momentum k when an emitting 
electron acquires the momentum transfer q. This probability has a universal 
form which is independent of properties of a target. For an electron traversing 
an amorphous medium this fact is reflected in formula (|5.7| ). Indeed, passing on 
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to a momentum space we have 
aduj 



dwr{q,, k) 
aduj 



dg riKQ{g) + r2K^{g) (1 - exp(-iq^)) 



TTU! 



+ r2F2 



In (x + v^TT 

1 ^ ; 



X 



Foix) 



2x' 



x^ 



X 



L In fx + Vl + xA - 1. 



x^ 



(5.17) 

Remind that q is measured in electron mass. The probabihty of radiation in 
this form was found in For a differential probability of scattering (here we 
consider the multiple scattering) there is a known formula (cp ( |2.5| ), ( ^.6| ) and 

(H) 



dws{q) = Fs{q)d'^q, F^(q) = j^-^ J d'^gexp (-iq^) exp {-Vs{g)l), 



Vs{g) = n J d'^q{l- exp{-iqQ)) a{q), 



(5.18) 



where c(q) is the cross section of single scattering. 



Using the formula ( |5.17| ) one can easily obtain to within logarithmic accuracy 
expressions ( ^.15| ),( ^.14| ). Both a radiation of boundary photons and a radiation 
in a thin target may be considered as a radiation of coUinear photons (see e.g. 
[ ^ ) in the case when an emitting particle deviates at large angle {{}s ^ 1/7, 9 3> 
1). Using ( ^.171) at X ^ 1 we find 



dwr{q) 



aduo 



d qdWr{q)Fs{q) 



ri + r2 (in q^ 

aduj 



T^UJ 



ri + r2 



(in; 



(5.19) 



For a thin target value of q^ is defined by mean square of multiple scattering 
angle on a thickness of a target /, and for boundary photons is the same but 
on the formation length If. However, if we one intends to perform computation 
beyond a logarithmic accuracy, the method given in this Section has advantage 
since there is no necessity to calculate -Fs(q) and in our approach a problem of 
calculation of the Coulomb corrections is solved in a rather simple way. 

6 A target of an intermediate thickness / ^ Ic 

It appears that used in Section 4 approach permits one to consider an important 
case when / ~ /c- 
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According to the partition of integrals over time in formula ( ^.4] ) into four 
domains we can write the probabihty of radiation as 

— = V ^ ^ = -Re + r,m (6 1) 

The integrals in I^'"^^ we compute on the assumption: z/q ^ 1, T ^ 1, z/qT ~ 
1, K = 1. Since integrals in I^^^ don't contain the logarithmic divergence, only 
the domain 4 contributes. In the domains 1-3 one of the integrals in J^^^ contains 
an integration over an interval < t < T and due to this reason dwi^2,i, oc T ^ 1. 
So, we consider /j^-* 



oo 

(1) 



j dti J dt2exp{-i{ti + t2)){0\exp{-i{Hoit2 - T))) exp{-iHT) 

— oo T 

oo oo 

X exp{iHoti) - exp{iHo{t2 - ti))|0) = J dh J dta exp(-i(ti + ta + T)) 



X J d'g, J d'g2Ko{0, Q„h) [K,{q„ Q^,T) - Ko{g„ Q^,T)]Ko{Q^,0,t2). 

(6.2) 

Here a calculation of integrals over and Q2 may be performed e.g. in a such 
way: 

• an integral over relative angle between and ^2 gives Jq{[3qiQ2) where 

Jnix) is the Bessel function, B = Bc = ; — — and B = Bn = — - for the 

^ ' ^ H He 2sinhz/T 2T 

first and second terms in the square brackets in the right-hand side of ( |6.2|) , 

• the remaining integrals over qi and Q2 can be found in tables. 

So, we have 

00 00 ^ 



4^^ = ^J dti j dt2 exp {-i{ti +t2 + T)) 




N{t,,t2 



h + t2 + T 



(6.3) 



(1 + z^^ti^a) sinh uT + z/(ti + ^2) cosh uT 



For z/q ^ 1 the contribution into integral the term with A^(ti, ^2) is of the order 
of I/z/q and this term may be neglected. With allowance for T <^ 1 we find 



00 00 



1 /" /" 1 

/f ^ = J j dt2exp (-2(ti + ^2)) ^ 



00 ^ + ^2 

00 X 1 °° 1 (6-4) 

f^e- fdy = -^ fe- = ^, 
8m J X J Am J An 
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where x = ti + t2, y = ti — t2. 

The contribution of the domain 4 into the term with r2 {if^ in ( |6.1|) ) contain 
two additional operators p (see (|4.4| )) which resuh additional factor ^1^2 ; 



m 



the integrand. Integration over the relative angle between and Q2 gives here 
Ji{I3qiQ2) and subsequent evaluation of integrals is similar to those for ( |6.3|) . We 
find 



00 00 



r(2) 



~ J dti J dt2 exp +t2 + T)) 







{ti + t2 + ry 



(6.5) 

The contribution into the integral with N'^{ti,t2) gives a domain 

^1,^2 ~ l/'^o ^ 1- Since T -C 1 as well, we can put an exponent in this integral 

equal to 1. So the integral is 



00 00 



00 00 



j dti J dt2N^{ti,t2) = J dx J dy 











{axy + b{x + y) + ay 



= 21n- = 21ncothz/T, 



(6.6) 



where x = uti, y = i't2, a = sinhz/T, b = coshz/T. The second term in (|6.5|) is 
calculated as (see (|6l^)) 



00 00 ^00 
J dti J dt2 exp {-i{ti +t2 + T)) ~ J dte 





00 



{t + ry 



dte~ 



t 



71 



(6.7) 



. - l~-C-lnT-l + z-, 

(t + T) 2 ' 



where t = ti + 12. Putting together ( |6.6|) and ( p.7|) we have 



iP = 21ntanhz/T - C - InT - 1 + i 



(6.8) 



For computation of /{ = 13 we will use Eq. (^^ 
/f ^ = ^Jdti I dt2 exp {-i{ti + ta)) ' ^ 







{uti + vt2)^ (sinh z/t2 + I'ti cosh i^t2)^ 

(6.9) 



Integrating by parts over ti with regard for exp(— 2^2) — 1 we have 

T 



r{2) 

II ' ^ 







■ 1 


-1 


dt2 




An J 




Vt2 









00 


T 




/*' 


/* 


2 exp 











1 



^{^1 + ^2) cosh z/t2(sinh i^t2 + ^^i cosh i^t2y 



(6.10) 
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The second of these integrals is proportional to T and can be neglected. The 
first integral gives 



(2) 



r(2) 



A-K tanh uT 



For a calculation of I2 we use formulae ([4 .51) and ( [4 .101) 

o T 



r(2) 
-'2 



47r 



j dt{t - T) 



1 



1 



sinh ft 



V 



dt 



coth vt — 



vt 



1 , sinh vT 

-r- In 7^ — ■ 

47r vT 



(6.11) 



(6.12) 



Combining all the contributions of four domains we obtain finally 



dw ^ 4a / ^^^1) ^ ^^^(2)\ ^ ^j^g r ^ ^ (ln(z/sinh i^T) - 1 - C) rs]. 



n=l 



(6.13) 

In the used units (T = a//2) the formation length ( |2.1| ) is (see also ( |2.5| ) and 

(iHD) 

qI = —T- (6-14) 



ale 

T 



1 



1^0 + 1 

In the case of thick target (T ^ tc, z/qT ^ 1) we have from ( |6.13| ) 

aal a e 



dw a r /I 
— ~ — [n + (In z/Q - 1 



C-ln2)r2] + r2^, 

TTu; v2 TTu; 



27rco' 



27r72 e' 
(6.15) 

This formula gives the probabihty of radiation at z/q ^ 1 (see (|2.28|) , ( p.43| )) 
where the contribution of boundary photons ( [4.14|) is included. 

In the case of thin target z/qT ^ 1 but when z/gT ^ 1 we have from ( |6.13|) 
the probability ( |5.15|) without term oc 1/ Lt. So we have {vqT = 4k) 



dw 
duj 



a 



1-C r2 



(6.16) 



Note, that when the value of the parameter uq is not very large, the accuracy 
of the formulae ( |6.15| ) and ( |6.16| ) may be insufficient. In this case one have to 
compute the next terms of the expansion, as it was done in Sections 4 and 5 
(see ( [4.14| ) and ( |5.15| )). The same is true for ( |6.13| ). A detailed analysis of the 
probability of radiation in the targets of an intermediate thickness will be carry 
out elsewhere. 
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7 A qualitative behavior of the spectral inten- 
sity of radiation 

We consider the spectral intensity of radiation for the energy of the initial elec- 
trons when the LPM suppression of the intensity of radiation takes place for 
relatively soft energies of photons: uj < uj^ e: 



= o 7 n In — , 



(7.1) 



see Eqs.( p.9| ), (|2.14|) , ( |2.15|) and ( p.36|) . This situation corresponds to the 
experimental conditions. 

A ratio of a thickness of a target and the formation length of radiation ( |2.1|) is 
an important characteristics of the process. If we take into account the multiple 
scattering and the polarization of a medium then the formation length ( |3.3| ) has 
the form 

27^ 



V CO 



-1 



(7.2) 



this ratio may be written as 



Piuj) = T (z/o + K 



+ 



Ijj. 



2 H 



T 



Tr = TiuJr 
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(7.3) 



rad 



where we put that ~ W — . Below we assume that uoc ^ which is true 

V uj 

under the experimental conditions. 



If P{uj, 



2Tc -C 1 then at u = tUc a target is thin and the Bethe-Heitler 



spectrum of radiation, which is valid at u ^ Uc 



dI{uo) 
duo 



constj will be also 

valid aX UJ < ojc in accordance with Eqs.(|5.7|) and ( |5.8| ) since Ak = z/qT = Tc <^ 1. 
This behavior of the spectral curve will continue with uo decrease until photon 
energies where a contribution of the transition radiation become essential. 

If P{uOc) ^ 1 (Tc ^ 1) then at uo > uOc a. target is thick and one has the LPM 
suppression for uo < uOc- There are two opportunities depending on the minimal 
value of the parameter /5. 



Pn 



-Tr 



UOr. 



UOr. 



P. 



2Tr 



UOr 



UOr 



2/3 



(7.4) 



If (3m <^ 1 then for photon energies uo > uoi it will be UO2 such that 



PM 



UOc 

UOo ~ 

^2 ^2 



(7.5) 
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and for uj < uj2 the thickness of a target becomes smaller than the formation 
length of radiation so that for lo <^lo2 the spectral distribution of the radiation 
intensity is described by formulae of Section 5. Under these conditions for 4A; = 
z/qT = Tc ^ 1 the spectral curve has a plateau 

dl 2aJ 

— = = const (7.6) 

duj 71 

in accordance with ( [5.13| ). Such behavior of the spectral curve (first discussed in 
[^ ) will continue until photon energies where one has to take into account the 
polarization of a medium and connected with it a contribution of the transition 
radiation. 

At (3m ^ 1 a target remains thick for all photon energies and radiation is 
described by formulae of Sections 2 and 3. In this case at -C cUc (t'o ^ 1) 
and uj ^ (up/uJc)^^^ Up (z/q ^ k) the spectral intensity of radiation formed inside 
a target is given by Exp. (|2.39|) and ( p.43|) and the contribution of the boundary 
photons is given by ( |4.14| ). 

Since a contribution into the spectral intensity of radiation from a passage 
of an electron inside of a target (oc T) is diminishing and a contribution of 
the boundary photons is increasing with u decrease, the spectral curve has a 
minimum at = Um- The value of Um may be estimated from equation 



hlnz/o + — =0, ^ IH 1= , 

24 z/o/ W2 4v^z/o 

4 UJ 





(7.7) 



When a value of Tc is high enough, the solution of Eq.( [7.7| ) doesn't satisfy the 
condition z/q ^ k and in this case the equation ( [7.7|) ceases to be valid. For 
determination of ujm in this case we use the behavior of the spectral intensity 
of radiation at k ^ z/q. In this case a contribution into radiation from inside 
passage of a target is described by ( |3.13D whilst the radiation of the boundary 




photons reduces to the transition radiation and its contribution is given by ( [4.11| ). 
Leaving the dominant terms {i^qT is a;- independent) we have 

+lnKj=0, ^ = 1, '^m = y, Um^d—i^p- (7.8) 

Since the value 7r^/12 ^ 0.8 is of the order of unity, the solution of ( [7^71) at 
Km ^ differs only slightly from ojm- Because of this, if the condition 
2Tc{ijJp/uJcY^^ ^ 1 is fulfilled, the position of the minimum is defined by Eq.(f77D. 

8 Discussion and conclusions 

Now we consider the experimental data ||T^-[T^ from a point of view of the above 
analysis. It is shown that the mechanism of radiation depends strongly on the 
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thickness of a target. So, we start with an estimate of thickness of used targets 
in terms of the formation length of radiation. From Eq. (|7.3|) we have that 

27ll I 

Tc = — — > 20 at - — > 2%. 

Oi-^rad J-'rad 

The minimum value of the ratio of a thickness of a target to the formation length 
of radiation is given by Eq.dTJ) {(3m ^ 2Tc{ujp/uJcf/^). For defined value of 
this ratio is least of all for the heavy elements. Indeed, the value of ujp = c^oT 

depends weakly on nucleus charge Z {uq = 30 -i- 80 eV), while Uc = — - — oc Z^. 

OiLrad 

Furthermore, the ratio ujp/ujc decreases with energy increase. Thus, among all 
targets with thickness / > 2%Lrad the minimal value of (3rn is attained for the 
heavy elements (W, Au, U) at the initial electron energy e = 25 GeV. In this 
case one has uj^ — 250 MeV, cjp ~ 4 MeV, (3rn > 2.5. Since the parameter 
is energy independent and the ratio Up/uc oc 1/e, the minimal value /3m > 5 is 
attained at the initial electron energy e = 8 GeV for all targets with thickness 
/ > 2%Lrad which can be considered as thick targets. 

As an example of obtained results we calculated the spectrum of the intensity 
of radiation in the tungsten target with thickness / = 2%Lrad at the initial 
electron energy e = 8 GeV and e = 25 GeV. The characteristic parameters of 
the radiation process for this case are given in the Table. We calculated the main 
(Migdal) term (Eq.(|]28D), the correction term (Eqs.(|]33D,(|]40D) taking into 
account an influence of the polarization of a medium according to (Eq.( |3.1I| )), 
as well as Coulomb corrections entering the parameters z/q (Eq.( |2.10| )) and L{gc) 
(Eq.( p^.35D ). The contribution of an inelastic scattering of a projectile on atomic 
electrons (quite small for the heavy elements) is not included although this could 
be done using Eq. (|2.45|) . We calculated also the contribution of the boundary 



photons Eq. ( [4. 121) . Here in the soft part of the spectrum uj < oJdi'^d — 2 MeV for 
e = 25 GeV) the transition radiation term (|4.11| ) dominates in (|4.12| ), whilst in 



the harder part of the boundary photon spectrum u > the terms depending 
on both the multiple scattering and the polarization of a medium ( |4.9| ) and ( |4.10|) 
give the main contribution; for e = 8 GeV we have ujd ~ 700 KeV . It is seen 
that we have for the boundary photons spectrum a smooth curve which eliminate 
difficulties mentioned in 0]. All these results presented separately in Fig. 2 as 
well as their sum (curve 5). Note, that for energy e = 25 GeV in the region of 
the minimum of the spectral curve 5 where the ratio of the target thickness to 
the formation length is minimal {f3m — 2.7, see Table) it may be that the target 
is not thick enough to use the formulae for a thick target. For a comparison 
with experiment we extract some data from Fig. 7 of [0. The theoretical curve 
gives the spectral distribution of the intensity of radiation (in units 2a /t^) without 



adjusting parameters. Data from were recalculated according with procedure 
given in it. One can see that agreement between the experiment and theory is 
rather satisfactory but far from being perfect. However, one has to take into 
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account that the theory of LPM effect in all previous papers had the logarithmic 
accuracy and did not contain Coulomb corrections. These shortcomings did not 
permit to pass to the Bethe-Heitler formula with acceptable accuracy and led to 
some difficulties in data processing. Both these shortcomings are overcome in the 
present paper. So, in our opinion, it is quite desirable to handle the experimental 
data using the formulae of this paper. 

The measurements in were made also using gold target with thickness 
I = 0.7%Lrad- For this case one has Tc ~ 6, /3m(25) = 0.7, /3m(8) = 1.5, so 
we have here a target of an intermediate thickness (see Section 6). We want to 
stress once more that for estimation of an effective thickness one have to use 
the formation length with regard for the multiple scattering and the polarization 
of a medium (see ( p.3|) and (|7!^) ). A detailed calculation for this case will be 
published elsewhere. 
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A Appendix 



A potential V{q) with the Coulomb corrections 

It is well known, that for heavy elements the Coulomb correction to the cross 
section of bremsstrahlung of high energy particles (correction to the Born ap- 
proximation) is quite sizable, see e.g. Eq.(18.30) in [0. The Coulomb correction 
(order of one for heavy elements) is subtracted from the "large" logarithm and 
if an accuracy of calculation goes beyond logarithmic one, one has to take into 
account this correction. For tungsten {Z = 74), gold {Z = 79) and uranium 
{Z = 92) in the case of complete screening the relative Coulomb corrections 
to the standard Bethe-Heitler cross section are respectively -7.5% , -8.3% and 
-10.7%. 

We consider the problem using eikonal approximation (see e.g. Appendix E in 
[0]). An amplitude /(q) and a cross section of scattering in this approximation 
have the form: 

/(q) = 77- f ^ exp (-iqe) 5(e), S{q) = exp {-ixig)) - 1, 

(A.l) 

x{q)= J U{Q,z)dz, da{q) = \f{q)\'d\ 

— oo 

where {z, g) are the longitudinal and transverse coordinates respectively, U{q, z) 
is the potential. Repeating a derivation made in Section 2 (eqs. (|2.3|) - (|2.9|) ) but 
with the cross section ( |A.1| ) we find for the potential V{q) 

V{Q)=nJ {l-exp{tqQ))\f{q)\'d'q 

, (A.2) 
= n d^x (S(x)5*(x) - S(x + q)S*{x)). 

Since the potential V{q) was calculated above in the Born approximation, we 
can calculate here the difference of the potentials calculated in the Born approx- 
imation Vb{q) and in the eikonal approximation V{q) 

AVig) = Vb{q) - V{q) =nj d^x{ exp [ix{q + x) 

oo ^ 

x(x) = / dz exp ( ) = 2ZaKo ( — ) , r 

J r \ KagJ 

— oo 

where Kq{z) is the modified Bessel function. Because the eikonal phase enters 
([A.3|) only in the combination x(e + x) — x(x) in the interesting for us region x ~ 



- ix{x)] - 1 

(A.3) 

= vz^ + x^, 
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X 



Q. Since Kq — ) is large only if x/a^ ^ 1 it is evident that main contribution 

\asJ 

into integral (|A.3| ) gives the region x ^ q <^ a^. In this region one has 



X(e + x) - x{x) = ^In' 



(A.4) 



So, in the expression for ^V{q) enters only one dimensional parameter q = Ig, 
where 1 is the unit vector. After substitution of variables x — > we have 



x + 1) 



2\ 



x^ 



2 



Changing variables y = ^ and then z = y + 1 we have 

X'^ 



211 J (z-l)4 



— In z 



(A.5) 



(A.6) 



Integration over azimuthal angle gives 

2lT 







{z"^ - 2zcos0 + 1)^ 



27r(l + z^) 

b2-l|3 • 



(A.7) 



Changing the variable = m, splitting the integration interval into two parts: 
(0, 1) and (1, oo) changing in the second interval f = 1/m we obtain 



Integrating by parts and changing once more variable u = e ^ we find 



fiO = Re 



(A.9) 



Integrating once more by parts and using the standard (Gauss) representation 
of the Euler ^/^-function we have finally 



AV{g) = 27rg^ f{Za), /(O = ^'Re [.^(1 + zO - ^(1)] = E 



^1 n{n^ + C^) 

(A.IO) 

The obtained function /(^) is the known Coulomb correction to the Bethe-Heitler 
cross section of bremsstrahlung, see e.g. Sections 17,18. 
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B Appendix 

An allowance for a form factor of a nucleus 



When Qc -C Rn (see (|2.15| ), (|2.16|) ) one cannot consider the potential of 
a nucleus as a potential of a point charge. A contribution into the multiple 
scattering gives a momentum transfer q < 1/ Rn- Because of the same reason the 
phase in expression (|A.2|) for the potential V{q) is small qg < Qc/Rn ^ 1 
and one can expand the potential. As a result we obtain 

.2 ^ J2 



ng 



J \cif{q)\'d'q = ^J |V5(x)|Vx 



|/(Vx(x)fci^a; = |/(q(x))^d\, 



(B.l) 



where q(x) is the classical momentum transfer on a straight-line trajectory with 
an impact parameter x. As one can see from ( |B.1| ), the mean square of the 
momentum transfer is the same the in eikonal approximation, in the Born ap- 
proximation and in the classical theory. The Coulomb correction in this case 
vanishes. 

Considering a nucleus as an uniformly charged sphere with the radius Rn, we 
have 



q(^) 



2^Q 



Q 



§{Q-Rn)+^ 



Rl 



d{Rn~Q) 



(B.2) 



\/l — X + xVl — X. 



Substituting the expression obtained into ( p.l|) we find the potential V{q) under 
conditions considered 







-1 - 2C) - 
Vig) = vreW 



7 



' Rn 



41n2 



In 



2a, 

Rn 



In 



Rl 



- 0.0407 



ln-^ + ^-2(C + ln2) 



If one uses standard representation of nuclear form factor (see e.g. |r7[ ) 
1 9 Ri 



^^^^ (1 + gW 

then one obtains 

J q'^{g)(fg = ATxZ'^a 



9o = ^, 



1.2 ■ 10"^^^/^ 



cm, 



ln-^ + ln6-2 
Ri 



2^,2 



ln^-0.20^ 



(B.3) 
(B.4) 

(B.5) 
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Taking into account that In ~ 20 we see that the difference between different 
models of nucleus is less than 1 % . 
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Figure captions 



Fig.l The functions Di 2(1^0) (Eq. (|2.4CI| )) vs parameter uq. 

Fig. 2 The intensity of radiation lu— — in tungsten with thickness / = 

diu 

2a 

0.088 mm in units — , ((a) is for the initial electrons energy e = 25 GeV 

and (b) is for e = 8 GeV). The Coulomb corrections and the polarization 
of a medium are included. 

— Curve 1 is the contribution of the main term ( p.28| ); 

— curve 2 is the correction (|2.33| ), ( p.40| ); 

— curve 3 is the sum of the previous contributions; 

— curve 4 is the contribution of the boundary photons (f4.12|) ; 

— curve 5 is the total prediction for the intensity of radiation. 
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TABLE Characteristic parameters of the radiation process in 
tungsten with the thickness I — 2%Lrad 



e {GeV) 


{MeV) 


Up {MeV) 




ui {MeV) 




uj^ {MeV) 


25 


228 


3.93 


21.25 


1.6 


2.7 


2 


8 


23.35 


1.26 


21.25 


0.76 


5.7 


0.5 
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Fig. 1 




1 10 100 

Photon energy CL>(MeV) 



Fig. 2 (a) 
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Fig. 2 (b) 



